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Abstract 

m ■ 

Two motions of oscillation and vacillating breathing (swing) of a red blood cell have been observed 

in bounded Poiseuille flows (Phys. Rev. E 85, 16307 (2012)). To understand such motions, we have 

£S| ■ studied the oscillating motion of a neutrally buoyant rigid particle of the same shape in Poiseuille 

flow in a narrow channel and obtained that the crucial point is to have the particle interacting 

O | with Poiseuille flow with its mass center moving up and down in the channel central region. Since 

.^^ the mass center of the cell migrates toward the channel central region, its oscillating motion of the 

inclination angle is similar to the aforementioned motion as long as the cell keeps the shape of long 

body. But as the up-and-down oscillation of the cell mass center damps out, the oscillating motion 

of the inclination angle also damps out and the cell inclination angle approaches to a fixed angle. 

Keywords Oscillating motion, red blood cell, neutrally buoyant particle, bounded Poiseuille 
^»' flow, narrow channel. 
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5h ■ 1 Introduction 

The microhydrodynamics of dcformable entities such as lipid vesicles and red blood cells in flows has 
received increasing attention experimentally, theoretically, and numerically in recent years. Lipid vesicles 
(e.g., see Q~]-[TT]) and red blood cells (e.g., see [H]-[I7]) show phenomenologically similar behaviors in 
Oh| shear flows, which are (i) a tank treading rotation with a stationary shape and a finite inclination angle 

with respect to the flow direction, (ii) an unsteady tumbling motion, and (iii) vacillating breathing (the 
long axis undergoes oscillation about the flow, while the shape shows breathing). In [18 ] 119]. two motions 
of oscillation and vacillating breathing of a red blood cell in bounded Poiseuille flows have been observed 

C-"- ' in low Reynolds number regime. The vacillating breathing motion is actually a combination of oscillation 

and deformation of the cell. The motion of oscillation of a blood cell in bounded Poiseuille flows is more 
interesting. In |20j , the motion of an elliptical cylinder particle immersed in an incompressible Newtonian 

f—^ , fluid in a narrow channel has been examined numerically in Stokes flow regime under the assumption 

that no external forces or torques act on the elliptical cylinder. Two interesting motions of a particle of 
elliptic shape are that it either tumbles (i.e., keep rotating) while the mass center always stays away from 
the centerline or rotates changing itself direction as the mass center crosses the centerline (oscillating 
motion). Even though the oscillating motion of a particle of elliptic shape is a periodic motion in Stokes 
flow regime, it closely resembles to the one associated with a red blood cell in [T5HT5] . 

In this paper we have compared the motions of a neutrally buoyant particle and a red blood cell in 
Poiseuille flow in the low Reynolds number regime to understand the oscillating motions in [TSJ Qjl] and 
to find out the difference between the cell motion and the particle motion in a narrow channel. For the 
motion of a neutrally buoyant particle of either biconcave or elliptical shape, we have obtained similar 
oscillating motion at the channel central region in Poiseuille flow when the particle mass center is placed 
at the centerline initially; but such motion is a transition since, later on, the particle migrates away from 
the centerline and starts tumbling as expected. Indeed in bounded Poiseuille flows (e.g., see [21]- [27]), a 
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neutrally buoyant particle migrates laterally to its equilibrium position between the channel centerline 
and the wall due to the competition between the shear gradient of the Poiseuille flow and the wall effect. 
For a neutrally buoyant cylinder of elliptical shape, it tumbles while migrating toward and then staying at 
its equilibrium position (e.g., sec [27]). Concerning the motion of a red blood cell, after the cell migrates 
to the central region the oscillating motion occurs due to the interaction between the cell and the profile 
of the Poiseuille flow just like an elliptic shape particle in [20] if the cell can maintain a long body shape. 
Specifically during the oscillating motion for the both cases, both mass centers move up and down in 
the channel central region. But unlike the periodic up-and-down motion of the elliptical cylinder mass 
center in Stokes flow in |20j , the cell oscillating motion damps out while its mass center approaches to an 
equilibrium height at the central region in a narrow channel and the cell inclination angle approaches to 
a fixed angle. For the cases of a neutrally buoyant cylinder of either biconcave shape or elliptical shape 
in the low Reynolds number regime, its mass center moves up and down about the centerline with an 
increasing amplitude and then moves away from the central region and toward its equilibrium position 
between the channel centerline and the wall. Thus the up and down motion of the mass center in the 
channel central region triggers the oscillation motion of a long body entity in Poiseuille flows. 

The scheme of this paper is as follows: We discuss the models and numerical methods briefly in 
Section 2. In Section 3, we have investigated the effects of the channel height, the bending rigidity of 
the membrane, the initial position on the motion of a red blood cell and a neutrally buoyant particle in 
Poiseuille flow in a narrow channel. The conclusions are summarized in Section 4. 

2 Models and methods 

2.1 Model and method for a red blood cell 
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Figure I: Schematic of a single red blood cell in plane Poiseuille flow with the computational domain Q. 

An RBC with the viscosity of the cytoplasm same as that of the blood plasma is suspended in a fluid 
domain fi filled with blood plasma which is incompressible and Newtonian as in Figure [T] For some 
T > 0, the governing equations for the fluid-cell system are the Navier-Stokes equations 

p f (-^ + u- Vuj = -\7p + pAu + f, infix (0,T), (I) 

V-u = 0, infi x (0,T). (2) 

Equations ([!} and (0) are completed by the following boundary and initial conditions: 

u = on the top and bottom of 51 and u is periodic in the x direction, (3) 

u(x,0) =u (x), in (4) 

where u and p are the fluid velocity and pressure, respectively, pf is the fluid density, and p is the fluid 
viscosity, which is assumed to be constant for the entire fluid. In eq. ([I}, f is a body force which is the 
sum of f p and %b where f p is the pressure gradient pointing in the x direction and fs accounts for the 
force acting on the interface between fluid and cell. In eq. ([3]), Uq(x) is the initial fluid velocity. 




Figure 2: The clastic spring model of the RBC membrane. 

The deformability and the elasticity of the RBC are due to the skeleton architecture of the membrane. 
A two-dimensional elastic spring model used in |28j is considered in this paper to describe the dcformablc 
behavior of the RBCs. Based on this model, the RBC membrane can be viewed as membrane particles 
connecting with the neighboring membrane particles by springs, as shown in Figure [2] Energy stores 
in the spring due to the change of the length I of the spring with respect to its reference length ^o 
and the change in angle 9 between two neighboring springs. The total energy of the RBC membrane, 
E = Ei + Et,, is the sum of the total energy for stretch and compression and the total energy for the 
bending which, in particular, are 

h N i _ ; i N 

i=l ° i=l 

In equation ([5)1. N is the total number of the spring elements, and ki and kb are spring constants for 
changes in length and bending angle, respectively. The cell shape is stimulated by reducing the total 
area of the circle of radius i?o = 2.8 /jm through a penalty function 

r s = ^(___-)2 (6) 

2 s e 

where s and s e are the time dependent area of the RBC and the specified area of the RBC, respectively, 
and the total energy is modified as E + T s . Based on the principle of virtual work the force acting on 
the ith membrane particle now is 

where r 2 is the position of the ith membrane particle. The value of the swelling ratio of an RBC in this 
paper is defined by s* = s e /(i:R^). 

The motion of the RBCs in the fluid flow is simulated by combining the immersed boundary method 
[^HHSHIM] and the aforementioned elastic spring model for RBC membrane. The Navier-Stokes equations 
for fluid flow have been solved by using an operator splitting technique and finite element method [321 I33j 
with a regular triangular mesh so that the specialized fast solver, such as FISHPAK by Adams et al. |34j . 
can be used to solve the fluid flow. The methodologies have been validated in previous studies [HI [33] . 

2.2 Model and method for a neutrally buoyant particle 

We suppose that fl is filled with an incompressible viscous Newtonian fluid of the density pf and 
viscosity /i. Let B(t) be a freely moving rigid neutrally buoyant particle in a fluid as in Figure [31 The 
boundaries of Q, and B{t) are denoted by T and OB, respectively. For some T > 0, the governing 
equations for the fluid-particle system are the Navier-Stokes equations for the fluid flow: 
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Figure 3: Schematic of a neutrally buoyant particle in plane Poiseuille flow with the domain tt. 



,du 



p f (— + (u • V)u) = f - Vp + 2^V • D(u) in fl\B(t), t & (0, T), 
V ■ u = in to\B(t), t e (0, T), 



u(x, 0) = u (x), Vx e O\B(0), with V • u = 0, 

u = on r, 

u = V p + ljx Gx, Vx e dB(t), 



(8) 

(9) 
(10) 
(11) 
(12) 



where u is the flow velocity, p is the pressure, 2D(u) = Vu + (Vu)*, and f is the pressure gradient 
pointing in the x direction. In (|I2p . the no-slip condition on the boundary of the particle, V p is the 

translation velocity, wx Gx= (—U)(y — G2),oo{x — Gi)) f where u is the angular velocity, G = {Gx^G^f 
is the mass center and x = (x, y)* is a point on the boundary of the particle. 
The motion of the particles is modeled by Newton's laws: 
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dt "'° ' '' dt ' ' dt ' Pl dt 
G(0) = G°, V p (0) = V°, lj(0) = w°, 6{0) = 8°. 
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In (fT3)) - (ri4|i . g is the gravity, 9 is the inclination angle between the long axis and the horizontal direction, 
M and I are the the mass and the moment of inertia of the particle, respectively; F and F* denote, 
respectively, the hydrodynamic force and the related torque imposed on the particle by the fluid given 

by 

F = — / crnds, F* = — / Gx x crnds, 
Job Job 

where <r is the stress tensor, x is the generic point on the boundary of the particle, n is the unit normal 

vector on the boundary of the particle pointing to the center of the particle and a x b = a\bi — 0261 for 

the two dimensional cases. 

The method of solution for the above fluid-particle interaction is a combination of a distributed 

Lagrange multiplier based fictitious domain method, the operator splitting methods and finite element 

methods (e.g., [27l [32 [35]). The basic idea is to imagine that the fluid fills the entire space inside 

as well as outside the particle boundaries. The fluid flow problem is then posed on a larger domain 

(the "fictitious domain"). This larger domain is simpler, allowing a simple regular mesh to be used, 

which in turn renders use of specialized fast solution techniques as in the previous section. The larger 

domain is also time independent, so the same mesh can be used for the entire simulation, eliminating 

the need for repeated remeshing and projection. The fluid inside the particle boundary must exhibit 

a rigid body motion. This constraint is enforced by using the distributed Lagrange multiplier, which 

represents the additional body force per unit volume needed to maintain the rigid body motion inside 

the particle boundary, much like the pressure in incompressible fluid flow, whose gradient is the force 

required to maintain the constraint of incompressibility The computational method has been validated 

in our previous studies j35j and |27j for the motion of neutrally buoyant disks and elliptical cylinders, 

respectively, in Poiseuille flow. In this paper, we have applied the method to the cases of neutrally 

buoyant cylinder of a long body shape freely moving in bounded Poiseuille flows. 



3 Simulation results and discussions 

3.1 Oscillating motion of a single RBC in a narrower channel 

Two motions of oscillation and vacillating-breathing (swing) of an RBC are observed in a narrow (100 
x 10 /im 2 ) channel considered here. The values of parameters for modeling cells are same with p"8l[T9] 
as follows: The bending constant is kb — 5 x 10~ 10 Nm, the spring constant is fc; = 5 x 10 -8 Nm, and 
the penalty coefficient is k s = 10 -5 Nm. The swelling ratios of the cells in the simulations are s* = 0.481 
and 0.9. The cells are suspended in blood plasma which has a density p = 1.00 g/cm 3 and a dynamical 
viscosity /i = 0.012 g/(cms). The viscosity ratio which describes the viscosity contrast of the inner 
and outer fluid of the RBC membrane is fixed at 1.0. The computational domain is a two dimensional 
horizontal channel. In addition, periodic conditions are imposed at the left and right boundary of the 
domain. The Reynolds number is defined by Re = pUH/p, where U is the average velocity in the 
channel, and H is the height of the channel. To obtain a Poiscuillc flow, a constant pressure gradient 
is prescribed as a body force so that the Reynolds number of the Poiseuille flow without cell is about 
0.4167. We consider the motion of a single RBC with four different bending constants which are O.lkb, 
lkb, lOfcf,, and 100kb in a Poiseuille flow with the fluid domain 100 x 10 /mi 2 (the degree of confinement 
Ro/w = 0.56). The initial velocity is zero everywhere. The grid resolution for the computational domain 
is 64 grid points per 10 /mi. 

First we have studied the cases in which the initial position of the cell mass center is located at the 
centerline and the initial inclination angle is 6 = it/ A. Different motions led by the different bending 
constants are observed. When the bending constants are O.lfcb and lkb, the cells of both swelling ratios 
stay at the center region of the channel and the parachute shapes have been obtained for both cells as 
shown in Figures [4] (a)- (d). For the bending constant 10kb, the cell of swelling ratio s* = 0.481 exhibits 
a damped oscillation with deformation (called the vacillating-breathing behavior [4] ) until it attains the 
equilibrium state aligning itself at an angle with the direction of the flow as shown in Figures 0] (e) and 
[51 The cell has a slipper shape whose mass center is slightly away from the center line as studied in [19] . 
But the one of swelling ratio s* = 0.9 only has oscillating motion for shorter period of time (see Figure 
S](/)) and the cell gradually deforms into a parachute shape. When the bending constant is 100fc&, both 
cells first exhibit damped oscillation with the shapes of long body and, at the end, align themselves with 
a fixed inclination angle with respect to the flow direction as in Figures 0] (g), (h) and [5] The histories 
of the height of the cell mass center in Figure [5] do show the correlation with the oscillating motion. The 
initial inclination angle 9 = 7r/4 helps the fluid flow to create the oscillation of the cell mass center and 
then the flow field, which is about a full quadratic profile, interacting with the long body creates the 
oscillation similar to the one in j20j . But the cell mass center gradually approaches to an equilibrium 
height due to the interaction between the deformability of the cell and the Poiseuille flow and the wall 
effect as in Figure [5] and at the same time the oscillating motion is damping out accordingly as shown 
by the histories of the inclination angle in Figure [5] Once the oscillation of the cell mass center damps 
out, the inclination angle then is fixed as in Figures \S\ and |5J 

When increasing the channel height to 20 /im and keeping the other parameters same, we have focused 
now on the cases of the bending constant lOOfcj,. In Figure [7j the histories show that the oscillation of 
inclination angle become periodic. These cells with the bending constant 100kb stay at the channel 
central region since they are not neutrally buoyant rigid particle and and the deformability plus the 
Poiseuille flow profile keeps them staying in the channel central region. But the wall effect is weaker in 
a twice wider channel so that the up and down motion of the mass center does not damps out fast as in 
the channel of height 10 /im. 

When placing the initial position of the cell mass center off the centerline in a channel of height of 10 
/mi , we have obtained almost similar behaviors for both cells in Poiscuillc flow as shown in Figures [51 [5] 
andlTUl The vacillating-breathing behavior for the case of the bending constant lOfcj, and the oscillating 
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Figure 4: (Color online). The snapshots of the cell migration in Poiscuillc flows for ,s*=0.481 and 0.9 
with different bending constants at different time (ms): (a) s* = 0.481 and O.lfcb, (b) s* = 0.9 and O.lkb, 
(c) s* = 0.481 and lk b , (d) s* = 0.9 and lk b , (e) s* = 0.481 and 10fc&, (f) s* = 0.9 and 10k b , (g) s* = 
0.481 and 100fc b , and (h) s* = 0.9 and lOOfcf,. The red asterisk denotes the same node point on the cell 
membrane. Rq/w = 0.56. The initial position is (5,5) at the channel centerline and the initial inclination 
angle is 7r/4. 
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Figure 5: (Color online). The history of the position of the cell mass center in Poiseuille flows for 
s*=0.481 and 0.9 with different bending constants: 10fc b (left) and 100fc b (right). R /w = 0.56. The 
initial position is (5,5) at the channel centerline and the initial inclination angle is tt/4. 
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Figure 6: (Color online). The history of the inclination angle 8 between the long axis of cell and the 
horizontal line for the bending constants lOkb (left) and lOOkb (right). The initial position is (5,5) at 
the channel centerline and the initial inclination angle is 7r/4. 
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Figure 7: (Color online). The history of the position of the cell mass center (left) and inclination angle 
(right) for s*=0.481 and 0.9 with the bending constant lOOfcf, in a channel of height 20 /iin. The initial 
position is (5,10) at the channel centerline and the initial inclination angle is 7r/4. 
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Figure 8: (Color online). The snapshots of the cell migration in Poiscuillc flows for ,s*=0.481 and 0.9 
with different bending constants at different time (ms): (a) s* = 0.481 and O.lfcb, (b) s* = 0.9 and O.lfcb, 
(c) s* = 0.481 and lk b , (d) s* = 0.9 and lk b , (c) s* = 0.481 and 10k b , (f) s* = 0.9 and 10*%, (g) s* = 
0.481 and lOOfc^, and (h) s* = 0.9 and lOOfcb. The red asterisk denotes the same node point on the cell 
membrane. The initial position is (5,3) and the initial inclination angle is 0. 
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Figure 9: (Color online). The history of the position of the cell mass center in Poiscuillc flows for 
s*=0.481 and 0.9 with different bending constants: 10/cj, (bottom left), and lOOfcf, (bottom right). The 
initial position is (5,3) and the initial inclination angle is 0. 
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Figure fO: (Color online). The history of the angle between the long axis of cell and the horizontal 
line for the bending constants lOkb (left) and lOOkb (right). The initial position is (5,3) and the initial 
inclination angle is 0. 

motion for the case of the bending constant lOOkj, are stronger for the cell of s* = 0.481 since the lateral 
migration toward the channel central region enhances the oscillation. But another interesting result is 
the motion of the cell of s* = 0.9 for the bending constants 10fc{, and 100fcb as shown in Figures [5] (/) 
and (h). A similar motion called snaking motion in Poiseuille flow has been studied in the Stokes regime 
in [57]. 

Concerning the oscillating motion and vacillating-breathing motion in Poiseuille flow, the bending 
constant needs to be large enough with respect to the maximum velocity of the fluid flow u max so that 
the cell can not be deformed into a symmetric parachute in the channel central region, but maintain a 
long body shape to interact with Poiseuille flow. 



3.2 Oscillating motion of a neutrally buoyant particle in a narrower channel 

The motion of a neutrally buoyant particle of either biconcave or elliptical shape in a narrow channel 
(100 x 10 /mi 2 ) has been studied in this section. The biconcave and elliptical shapes are obtained based 
on the cell resting shapes used in the previous section. The long axis of the biconcave shape (resp., 
elliptical shape) is 7.65 pra (resp., 6.825/^m). The particle is suspended in an incompressible Newtonian 
fluid of the density p = 1.00 g/cm 3 and the dynamical viscosity /i = 0.012 g/(cms). A constant pressure 
gradient is chosen as in the previous section so that the Re of the flow without particle is about 0.4167. 
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Figure 11: (Color online). History of particle mass center and inclination angle of the biconcave particle 
with s* = 0.481 (top two) and the elliptic particle with s* = 0.9 (bottom two). 
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Figure 12: (Color online). Snapshots of particle motion in a narrow channel with the initial position of 
the mass center at the channel centerline. (i). Biconcave shape: 0° < 9 < 90° for 1.5 < t < 1.9 and 
710.7 < t < 713.7, 0° < 6 < 180° for 719.7 < t < 725.7 and 1144.1 < t < 1150.4, and tumbling for 
1499.2 < t < 1501.7. (ii) Elliptical shape: 0° < 6 < 180° for < t < 2 and 272.5 < t < 278 and tumbling 
406.3 < t < 407.6 and 957.2 < t < 958.1. The number above each particle indicates the time (time unit 
is ms). 




Figure 13: (Color online). The histories of particle mass position (top) and the angular velocity (bottom) 
with s* = 0.481 (left two) and s* = 0.9 (right two) for < t < 20. 
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Figure 14: (Color online). History of particle mass center and inclination angle of the biconcave particle 
with s* = 0.481(top two) and the elliptic particle with s* = 0.9 (bottom two) in a channel of height 20 
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Figure 15: (Color online). History of particle mass center of the biconcave particle with s* = 0.481 (left) 
and the elliptic particle with s* = 0.9 (right). The initial mass center is off the center line. 
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Figure 16: (Color online). Snapshots of particle motion in narrow channel with the initial position of 
the mass center away from the channel centerline. The number above each snapshot denotes when the 
snapshot was taken(time unit is ms). 
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Figure 17: (Color online). History of particle inclination angle when forced to stay at channel center 
with s* = 0.481 (left) and s* = 0.9 (right) 
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We have first considered the motions of the particle mass center of both shapes initially located at 
the centerline of the channel with the initial inclination angle 9 = ir/4. As in Figures [TT] and [T^l both 
particles show an oscillating motion of the mass center about the channel center. The up and down 
motion together with the Poiseuille flow velocity profile cause the particles to oscillate in the inclination 
angle. As the up and down motion becomes stronger, the oscillation in inclination angle keeps increasing 
its range. For the particle of biconcave shape, the oscillation of inclination angle is first between and 
90 degrees and then the particle starts to swing back and forth with the inclination angle oscillating 
between and i80 degrees. Once the particle turns horizontal, it starts a tumbling motion and ends 
the oscillation of the inclination angle. The mass center of the tumbling particle migrates toward an 
equilibrium position between the channel centerline and the wall. The initial inclination angle 8 = 7r/4 
helps the fluid flow to create the oscillation of the particle mass center as shown in Figure Q2J At the 
beginning once the mass center is pushing down by the Poiseuille flow, the inclination angle is decreasing 
since the upper portion of the particle close to the centerline moves forward faster than the lower portion 
away from the centerline does. The wall effect from the bottom wall slows down the lateral migration 
and finally push the neutrally buoyant particle away. When it crosses the centerline, the angular velocity 
changes sign accordingly as in Figure 1131 The neutrally buoyant particle goes through the same up 
and down motion with rotation changing its direction accordingly which is similar to the one in the 
Stokes regime in |20j . But the neutrally buoyant particle moves up and down about the centerline with 
increasing amplitude and finally breaks away from the centerline and moves to its equilibrium position 
between the centerline and the wall as shown in Figure 1111 The elliptic shape particle has a similar 
motion, but the oscillation of the inclination angle is always a back and forth swing motion before the 
tumbling motion occurs. This difference in motion is due to the shape difference of the particles. The 
particle of biconcave shape undertakes stronger resistance force from the flow while the flow can smoothly 
bypass the particle of elliptic shape. This also explains why it takes the particle of biconcave shape longer 
time to turn to a tumbling motion than the particle of elliptic shape does. 

When increasing the channel height to 20 /xm and keeping the other parameters same, both the up 
and down motion of the mass center and the rotation of the long body oscillates stronger since the wall 
effect is weaker in a twice wider channel. The neutrally buoyant particle of long body shape behaves 
similarly but they migrate away from the centerline faster than they do in the channel of height 10 /im. 

When the initial mass center is off the centerline, the particle directly turns into a tumbling motion 
and migrates toward the equilibrium position between the channel centerline and the wall in Figure 
1151 When the particle reaches its equilibrium position, it shows a rotation with periodically varying 
angular velocity as in |27| . The neutrally buoyant particle behavior is entirely different from the cell 
motion discussed in the previous section due to the lack of the deformability. In Figure HH two plots 
of snapshots of the tumbling motion are shown for each swelling ratio. The first plot is for the first 
tumbling motion of the particle where as the second plot is the snapshot of one tumbling motion when 
the particle has reached its equilibrium position. By comparing the time in the first and second plots, 
one can observe that it takes longer time for the particle to spin a whole circle at the equilibrium position 
than in the first few spins. This is because under the parabolic velocity profile, the flow shear rate is 
higher when the particle is further away from the channel center. So the angular velocity is larger in the 
first few spins since the position of the particle is closer to the channel wall. 

To further study the effect of the up and down motion about the centerline of the mass center, we 
have added a constraint on the particle motion so that the particle is only allowed to move freely in the 
horizontal direction and rotate freely. In Figure [T71 the inclination angle of the constrained motion of 
a neutrally buoyant particle reaches to an equilibrium angle without any oscillation when moving along 
the centerline in the channel. It shows that the up and down motion of the mass center in the central 
region of the Poiseuille flow causes the oscillation of the inclination angle. 
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4 Conclusion 

We have compared the oscillating motions of a neutrally buoyant particle and a red blood cell in 
Poiseuille flow in a narrow channel to understand the oscillating motions in [181 119) and to find out the 
difference between the cell motion and the particle motion. For the motion of a neutrally buoyant particle 
of either biconcave or elliptical shape, we have obtained oscillating motion in Poiseuille flow when the 
particle mass center is placed at the centerline initially. But the neutrally buoyant particle moves up and 
down about the centerline with the amplitude which is increasing in time and finally breaks away from 
the centerline and moves to its equilibrium position between the centerline and the wall. The neutrally 
buoyant particle behavior is entirely different from the cell motion due to the lack of the deformability. 
Concerning the cell motion in the central region of the channel, the oscillating motion occurs exactly like 
a long rigid body as long as the cell mass center moves up and down in the channel central region and 
the cell can maintain a long body shape. But when the mass center of the neutrally buoyant particle of a 
long body shape is not allowed to move up and down, its inclination angle reaches a fixed angle without 
any oscillation. Thus the up-and-down motion of the mass center in the channel central region triggers 
the oscillation motion of a long body entity in Poiseuille flows. 
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